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Abstract In large-eddy simulations of atmospheric boundary layer tur- 
bulence, the lumped coefficient in the eddy-diffusion subgrid-scale (SGS) 
model is known to depend on scale for the case of inert scalars. This scale de- 
pendence is predominant near the surface. In this paper, a scale-dependent 
dynamic SGS model for the turbulent transport of reacting scalars is im- 
plemented in large-eddy simulations of a neutral boundary layer. Since the 
model coefficient is computed dynamically from the dynamics of the resolved 
scales, the simulations are free from any parameter tuning. A set of chemical 
cases representative of various turbulent reacting flow regimes is examined. 
The reactants are involved in a first-order reaction and are injected in the 
atmospheric boundary layer with a constant and uniform surface flux. Em- 
phasis is placed on studying the combined effects of resolution and chemical 
regime on the performance of the SGS model. Simulations with the scale- 
dependent dynamic model yield the expected trends of the coefficients as 
function of resolution, position in the flow and chemical regime, leading to 
resolution-independent turbulent reactant fluxes. 
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1 Introduction 

The atmospheric boundary layer (ABL) is the lower part of the atmosphere, 
where turbulence plays an essential role on the transport of momentum and 
scalars (heat, water vapor and chemicals) between the biosphere and the 
upper atmosphere. Quantifying the transport rate of chemicals from the 
Earth's surface into the upper parts of the atmosphere is critical for a wide 
range of environmental applications, including air quality modeling and re- 
gional and global climate studies. One important open question is how to 
represent the effects of boundary layer physical processes, in particular tur- 
bulence, on chemical transport and transformation of chemical compounds. 
A large number of reactive gases emitted by vegetation, such as biogenic 
volatile organic compounds, have typical lifetimes that are of the same order 
of magnitude or smaller than the mixing time scales of the largest ABL eddy 
motions J^. For such species (isoprene and monoterpenes for instance) the 
interplay between turbulence and chemistry is of great importance. Accurate 
modeling of their chemical transformations in the ABL requires improved 
understanding on how turbulence and spatial distribution of the reactants 
affect chemistry. 

Realistic numerical experiments on the turbulent transport and mixing 
of reactants require the use of Large-Eddy Simulation (LES). In the atmo- 
spheric boundary layer, the scales associated with turbulent motions range 
from the Kolmogorov dissipation scale (on the order of a millimeter) to 
the boundary layer depth (on the order of a kilometer). The largest eddies 
are responsible for the turbulent transport of the scalars and momentum 
whereas the smallest ones are mainly dissipative. LES consists of explic- 
itly resolving all scales larger than the grid scale (on the order of tens of 
meters in the ABL), while the smallest (less energetic) scales are param- 
eterized using a subgrid-scale (SGS) model. Since the pioneering paper of 
Schumann that showed for the first time that segregation of reactants 
due the inefficient mixing of the convective turbulence plays an important 
role on moderating the reaction rates, several numerical studies on the effect 
of turbulence on chemistry in the ABL have been performed using LES [HMl 
|^]| ^ |[7||^l F?lll()lllllll2lli;^lll4lll5lll()lll7j . In particular, Gao et al. and Patton 
et al. J5 focused their study on the lowest part of the ABL and on the 
interplay between chemistry and turbulence in and above canopies. 

An important challenge in large-eddy simulations of the atmospheric 
boundary layer is the specification of the SGS model coefhcients and, in 
particular, how to account for external or internal parameters such as posi- 
tion in the flow, grid/filter scale, atmospheric stability and chemical regime. 
Usually in LES of atmospheric turbulent reacting flows, the subgrid reactant 
flux is modelled as the product of a local gradient of the concentrations from 
the resolved scales and of eddy-diffusity coefficients. These coefficients are 
either simply prescribed (in the case of Smagorinsky models) or estimated 
by solving an equation for the turbulent kinetic energy assuming that sub- 
grid scale turbulence is equally efficient at transporting scalars as it is for 
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heat JH|- In this paper, we exphcitly calculate the model coefficients by us- 
ing a dynamic procedure and without making any assumption of the ability 
of subgrid turbulence to transport scalars. Thus the model coefficients are 
specifically calculated for each reacting scalar. 

When studying the effects of the unresolved scales on the chemical trans- 
formations, e.g. related to the subgrid scale mixing between the reactants 
due to subgrid scale turbulence, Vinuesa and Porte- Agel ^\ suggest that 
these effects can be correlated with the chemical regime. Since the SGS 
reactant flux includes a chemical contribution in its governing equation, 
chemical regime dependency of the subgrid scale transport can be also ex- 
pected. By using a set of chemical setups representative of several chemical 
regimes, one can study the relevance of this chemical dependency of the 
reactant subgrid fluxes. 

Recently, Porte- Agel [201 showed that the lumped coefficient used in the 
eddy-diffusion models for inert scalars are strongly scale dependent, in par- 
ticular close to the surface. In order to account for that scale dependence of 
the coefficient, he developed a scale-dependent dynamic procedure to com- 
pute the value of the eddy-diffusion coefficient using information contained 
in the resolved scales. In this paper, we use this tuning-free scale-dependent 
dynamic model for simulations of reacting scalars. 

The structure of this paper is as follows. In section|21 we present a scale- 
dependent dynamic model for the turbulent transport of reactants in the 
ABL. The model is used for simulations of a neutral atmospheric boundary 
layer. The characteristics of the simulations are presented in section |3| A 
set of simulations is performed in order to study the combined effects of 
chemical regime and resolution. The relative importance of subgrid-scale 
turbulence on the chemistry is classified using a subgrid-scale Damkohler 
number. By focusing our analysis on a first order decaying scalar, we are 
able to isolate the effects of the unresolved scales on the turbulent trans- 
port of the reactants. Emphasis is placed on the scale- independence of the 
simulated reactant fluxes and on the combined effects of chemical regime 
and resolution on the dynamically computed model coefficients and scale- 
dependence parameters. Finally, conclusions are drawn in sectional 



2 The scale-dependent dynamic model 

The filtered governing equation for the concentration of a reacting scalar A 
that is solved in LES reads: 

where ca and are the spatially filtered (at the filter scale A) and the 
subgrid concentrations of the reactant, respectively. u"c^ is the subgrid- 
scale flux and Rch is the chemical term. For a reactant A involved in the 
first-order reaction 
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A 



Products, 



(2) 



the chemical term reads 



Rch = ~jCA- 



(3) 



The effect of the umesolved scales on the evolution of the filtered scalar 
concentration appears through the subgrid scale reactant flux 



It is important to note that in the case of a second-order reaction (not 
considered here), the unresolved scales can affect the reactant concentrations 
not only through the subgrid-scale fluxes but also through the subgrid-scale 
covariances (see e.g. 

The subgrid reactant flux is usually modeled as a function of a local 
gradient of the resolved concentrations using an eddy-diffusity model. The 
subgrid-scale Schmidt number in that model is either simply prescribed as 
a constant value or estimated assuming that subgrid scale turbulence is 
equally efficient at transporting reacting scalars as it is for heat. 

The eddy-diffusion subgrid model is widely used in LES of the atmo- 
spheric boundary layer. A common formulation of this model for the subgrid 
scale flux of a scalar -0 is 



solved strain rate tensor, A is the filter scale and C0 is the resolved (spatially 
filtered) scalar concentration. Cs and Scsgs,ip denote the Smagorinsky co- 
efficient in the eddy-viscosity model and the subgrid-scale Schmidt number 
for the scalar ■0, respectively. 

Usually the subgrid-scale Schmidt number Scsgs,ip and, as a consequence, 
the lumped coefficient Sc~g^ are only determined (calculated or pre- 
scribed) for inert scalars, and those values are used for the calculation of 
the subgrid scale reactant fiuxes using Equation (Q. It is, therefore, assumed 
that chemistry does not affect the subgrid scale fluxes and, consequently, 
Scsgs,'4, = Scsgs, wherc Scggs is the subgrid-scale Schmidt number for an 
inert scalar. Scggs is often taken as a constant value {Scggs ~ 1/3), which is 
well-established for isotropic turbulence 22,23 . However, near the surface 
of ABL flows, the reacting scalar field at the smallest resolved scales can 
become anisotropic due to the combined effects of proximity to the surface, 
resolution (filter/grid scale), atmospheric stability and/or chemical regime. 
In the following, we propose an alternative procedure to calculate the value 
of the lumped eddy-diffusion model coefficient that is able to adjust to the 



Uj = UiCA - UiCA- 



(4) 
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anisotropy of the reacting scalar field associated with the aforementioned 
factors. In particular, the value of the lumped model parameter Sc~g^ ip^'s 
is determined based on the dynamics of the smallest resolved scales us- 
ing the scale-dependent dynamic model introduced for passive scalars by 
Porte- Agel |20| . This tuning- free procedure is summarized below. 

For the reactant subgrid scale flux m-'c^ (denoted qA,i hereafter to en- 
hance readability) , the dynamic procedure is based on the Germano identity 
for scalars [^11^ 



KA,i = Qaa - qA,i = "iC^ - UiCA, (6) 

where QA,i — UjCA — UiCA is the subgrid scale flux at a test-filter scale 
(typically taken as Zi = 2A) and KA,i is a resolved reactant flux that can be 
evaluated based on the resolved scales. Applying the eddy-diffusion model, 
QA,i is determined by 



iA,i 



Sc:,\aC's (^) 



a' 



S 



dcA 



(7) 



Substitution of Equations. l(SJl and {Tj) into © leads to the system 



Ka,i = Sc^g^j^CsXA,z: (8) 

where, for A — 2A, 



Xa. 



A' 



dcA 
dxi 



Sc 



sgs,A 



Cl{A) 



S 



(9) 



Note that the traditional dynamic model ^20^ assumes scale invariance 
of the model coefRcient at the filter and test filter scales, i.e., 



Sc:'aCI (A) ^ Sc-'Cl {2A) Sc^l^, 



(10) 



Minizing the error associated with the use of the eddy-diffusion model 
in Equation ^ over all three vector components results in 



Sc 



s,s,aCI{A) = 



(KA,iXA,i) 
{XAaXAa) ' 



(11) 



where () denotes an averaging operator (typically over horizontal planes) 

and XA,i contains the scale-dependent parameter (3 a — gj-i"^ p2 ■ Please 

note that /Jyi = 1 in the case of the traditional dynamic model. A dynamic 
procedure can be developed to compute the value for Pa using a second 
test-filtering operation at scale A > A. For the sake of simplicity, we take a 
second test- filter scale as Zi = 4 A, and denote variables filtered at scale 4 A 
by a carets This results in a second equation for ScJ^^ {A) that reads 
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where 



and 



UtCA - UiCA, 



(12) 



(13) 



S 



d^A ^, Sc;^\^CU4A) 



s 



^ ■ (14) 



By combining Equations. Illl and ll2l one obtains the foUowing equation 
from which 13a can be computed (more details on this procedure are given 
in Porte- Agel 

{KA^^XA^^) {X'AA^) " (^^l^.^A,.) = 0. (15) 

The value of (3a is extracted from Equation 1151 assuming that 



Sc-'aCI{A) Sc-J.Cl{2A) 



(16) 



which implies that ^cj^ (4Z\) /5c;i [A) = Pj 



Then Pa is used in Equation^Jto obtain the lumped coefficient Sc^^^ a^s- 



3 Reacting atmospheric flow simulations 

3. 1 Turbulent reacting flow classification 

Turbulent reacting flows can be classified by using two dimensionless num- 
bers pil2f)j . The first one, i.e. the turbulent integral Damkohler number, Dat, 
refers to the influence of the largest atmospheric boundary layer eddies on 
reacting scalars. In the context of Reynolds Averaged Navier Stokes (RANS) 
modeling, it gives a measure of the relative importance of the dynamical and 
chemical contributions to the reacting scalar mean governing equation. The 
second number, i.e. the Kolmogorov Damkohler number, Da^, quantifies 
the relative magnitude of the characteristic time scales of the smallest ABL 
eddies with respect to the chemical reaction time scale. These numbers are 
defined as 

Dat = -, (17) 



Dak = (18) 
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where Tt is the integral turbulent time scale, Tk is the Kolmogorov time scale 
and Tc is the chemical time scale. 

Using these Damkohler numbers, the following classification of chemical 
regime in turbulent reacting flows can be made: 

— Dat < 1, slow chemical regime: the reactant is uniformly mixed in the 
boundary layer and its governing equation is dominated by the dynam- 
ical terms. 

— Dak < 1 < Dat, moderate chemical regime: the largest scales affect the 
turbulent mixing whereas the chemistry is not limited by the small-scale 
turbulence. 

— Dak > 1, fast chemical regime: turbulence controls the chemistry at all 
scales. 

In the context of large-eddy simulation, it is also very important to assess 
the effects of the subgrid scales on chemistry. In order to achieve that, it is 
convenient to introduce the so-called subgrid Damkohler number or Dasgs 

m 

Da^gs = (19) 

Dasgs is a measure of the relative magnitude of the time scale associated 
with the smallest eddies resolved in LES with respect to the chemical time 
scale. That time scale can be computed as 



' sgs 



(20) 



1/2 

Esgs , 

where Vsgs and esgs are the subgrid-scale eddy viscosity and the subgrid- 
scale dissipation rate of the turbulent kinetic energy, respectively. The subgrid- 
scale dissipation rate is defined as the rate of transfer of kinetic energy be- 
tween the resolved and the subgrid scales |27II28| . By using Dasgs , turbulent 
reacting flows are classifled with respect to the effect of the unresolved scales 
on the chemistry. For Dasgs << Ij the flow is well-mixed at subgrid scales 
and, consequently, the unresolved scales do not influence the behavior of 
reacting scalars. For Dasgs ~ 0(1) and > 1, all the scales of turbulence, i.e. 
resolved and subgrid, affect the chemistry. These effects on the transport 
and mixing of reactants lead to a heterogeneous distribution of the reactants 
at subgrid scales. 



3.2 Numerical set-up 

The code used is a modifled version of 3-dimensional LES code described 
by Albertson and Parlange [211, Porte- Agel et al. 123, and Porte- Agel PH] . 
Briefly, the code uses a mixed pseudospectral finite-difference method and 
the subgrid-scale stresses and momentum and temperature fluxes are param- 
eterized with scale-dependent dynamic models. Test filtering for dynamic 
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models is done using two-dimensional spectral cutoff filters. A chemical 
solver has been implemented in the LES code ^H!- We simulate a neutral 
atmospheric boundary layer driven by a pressure gradient forcing and a 
surface roughness length for momentum of 0.1m is prescribed. The compu- 
tational domain is of size (L^, Ly, Lz) and it correspond to a grid prescribed 
with N X N X N points. The height of the domain is equal to 1000 meters 
and Lx = Ly = 2'kLz- Several resolutions are used with N equal to 32, 64, 
and 128. 

In order to restrict our study to the performance of the scale dependent 
subgrid-scale model for the reactant flux, the studied reacting flow consists 
of the first-order reaction ^ . 

We define three chemical setups based on this first-order irreversible 
reaction. The reactant so-called bottom-up reactant, is uniformly emitted 
at the surface with a flux of 0.1 pph m and no initial concentrations in 
the boundary layer. The reaction rate j is set to zero, 4.5 x 10~^ , and 
9 X 10""^ . The corresponding chemical set-up will be refer to as inert 
(/), slow {S) and fast (F) chemical cases. Notice that, in our simulations, 
Tt — hju^ and Tc = lead to turbulent Damkohler numbers Dat equals 
to 0, 1 and 20, respectively. The simulations cover a 1.5 hours period and 
the statistics presented here are obtained averaging the results over the last 
hour of simulation. 



3.3 Results 

In Fig.Q] the vertical profiles of the subgrid scale Damkohler number Dusgs 
are shown. For the slow (S) chemical regime, the chemical time scale (tc) 
is much larger than the characteristic subgrid-scale time scale {rggs), which 
determines a relatively small value of Daggs for all resolutions. On the con- 
trary, all the fast {F) chemical cases have Daggs ^ 0{1) and, as a result, 
the transport of the reacting scalar is affected by the subgrid scale chem- 
istry. The value of Daggs and, consequently, the effect of the subgrid-scale 
chemistry on the scalar fluxes, increases with decreasing resolution due to 
the fact that the characteristic time scale of the subgrid fluxes is larger for 
coarser grids. 

Fig. 13 shows the vertical proflles of the time-averaged total scalar fluxes 
obtained as the sum of the resolved and the subgrid-scale fluxes. Within the 
boundary layer, the profiles of inert scalar fluxes have a linear shape. As no- 
ticed by Porte- Agel |^ , since a constant mean surface flux is imposed and 
no viscous effects are considered, the averaged total turbulent flux decreases 
linearly with height. The reactive scalar flux profiles, however, deviate from 
this linear behavior. In fact, these deviations become more significant when 
the turbulent Damkohler number increases; they are larger for the F ex- 
periment than for the S one. The deviations increase with the reaction rate 
and, thus, with the turbulent Damkohler number due to the increase in the 
chemical contribution to the flux. This was shown also by Gao and We- 
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Fig. 1 Vertical profiles of the subgrid scale Damliohler numbers. These numbers 

have been obtained dynamically and averaged over the last hour of simulation. 
The S and F chemical cases are represented by dashed and dash-dotted lines, 
respectively, and the 32^, 64^ and 128^ grid-points simulations are plotted with 
squares, circles and no symbol, respectively. 
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Fig. 2 Vertical profiles of the total scalar fiuxes averaged over the last hour of 
simulation. Results are made dimensionless by using the surface fluxes. 
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Fig. 3 Vertical profiles of /3a as a function of the height made dimensionless using 
the boundary layer depth h (upper figure) and the filter size A (lower figure). The 
results have been averaged over the last hour of simulation. 

sely Sykes el al. and Vinuesa and Vila-Guerau de Arellano for 
convective boundary layers. 

As mentioned earlier, the magnitude of the subgrid-scale fluxes are af- 
fected by both the resolution and the chemical regime. The fact that the 
total turbulent fluxes (resolved plus subgrid contributions) for each chemical 
regime show almost identical results for the three resolutions under consid- 
eration highlights the robustness of the scale-dependent dynamic model. 
This subgrid model is able to adjust the magnitude of the subgrid-scale 



Subgrid-scale reactant flux model 



11 



fluxes to yield the same total turbulent fluxes. It is important to note that 
the resolution independence of the results is one of the most challenging 
tests for any subgrid-scale model, especially when the subgrid-scale effects 
are also affected by the chemical regime. 

The time-averaged values of the scale-dependence parameter j3A are 
shown in Fig. O as a function of the normalized heights zjh and z/A. 13a 
is smaller near the surface, indicating a strong scale dependence in that 
region, and it increases to larger values far from the surface. Note, however, 
that the scale-dependence coefficients obtained far from the surface are still 
slighly smaller than the value of 1, indicating a slight scale dependence and 
flow anisotropy. This behavior contrasts with the value of 1 obtained for the 
scale-dependent parameter for the momentum fluxes as reported by Porte- 
Agel 1201, and it could be attributed to the more anisotropic behavior of 
scalars (compared with the velocity field) in turbulent flows. From Fig. |3 
it is also clear that (3a depends on both resolution (scale A) and chemical 
regime (characterized by the subgrid Damkohler number). The scale depen- 
dence of (3a was already demonstrated by Porte- Agel ^50' for the case of an 
inert scalar, by showing that the values of (iA computed at different resolu- 
tions collapse when plotted against z/A. The fact that in our simulations 
the collapse of the curves is not perfect is due to the additional effect of 
chemistry on the subgrid-scale fluxes. This effect is further explored below. 

Fig. ^ show the vertical profiles of scale-dependent parameter (3a and 
the lumped coefficient Sc~gg j^C^ as a function of the normalized height 
z/A. In order to isolate the effect of the chemical regime from that of the 
resolution, the results from each resolution are presented in different pan- 
els. One can notice that for both (3a and Sc~g^ a^'s^ results from the 
inert and slow chemical cases are the same for the highest resolution. They 
are not affected by the chemistry in the / and S chemical cases, showing 
only a scale (resolution) dependence. The little effect of the chemistry on 
the subgrid ffuxes in the slow chemical regime is due to the fact that the 
characteristic time scale of the chemistry is much larger than the charac- 
teristic time scale of the subgrid scale fluxes, as illustrated by the very low 
values of the subgrid Damkohler number. When decreasing (coarsening) 
the resolution, however, chemistry has a relatively more important effect on 
the subgrid-scale fluxes, which translates also into an increase in the subgrid 
Damkohler number. As a result, the deviation from the / chemical case pro- 
files increases when coarsening the resolution. Chemistry has a clear effect 
on all F experiments with this effect increasing as the resolution coarsens 
and, consequently, Daggs increases. In summary, our results show that the 
scale-dependent dynamic model is able to caputurc the dependence of the 
subgrid coefficient and scale-dependent parameter on both resolution and 
chemical regime. 

The scalar variance spectra of the reactive scalar obtained from all the 
simulations at a height of z = 0.3/i are shown in Fig.jS] The spectra are cal- 
culated from spatial information using one-dimensional Fourier transforms 
that are then averaged in the spanwise direction and also in time. They 
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Fig. 4 Vertical profiles of (3a (left column) and lumped coefficient Sc^^^ ^Cs 
(right column) for the 32'' (upper row), 64'' (middle row) and 128^ (lower row) 
simulations. The results have been averaged over the last hour of simulation. 



are normalized with cr^z, where is the area below the spectral curve, and 
they are plotted against fciz, where fci is the streamwise wavenumber. Doing 
so allows us to focus on the distribution of the variance over the range of 
resolved scales. |2I3 showed that the normalized spectra for inert scalars and 
heights z ^ 0.3h collapse and are proportional to ki in the inertial sur- 
brange, which corresponds to scales smaller than the height (i.e., kiz > 1). 
At the largest scales (smaller wavenumbers), our results agree with Jonker et 
al. ^nii who showed that the effect of the chemistry on the variance spectra 
depends on the chemical regime, characterized by the turbulent Damkohler 
number Dat- For slow chemical regimes, where Dat << 1, the time scale of 
the chemistry (tc) is much slower than the integral time scale of the turbu- 
lence (rt) and, consequently, the chemical reaction has practically no effect 
on the scalar variance at any of the turbulence scales. When Dat ^ 0{1), 
the chemistry is faster than the characteristic turning time of the largest 
eddies, which has an impact (reduction) on the level of fluctuations of the 
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Fig. 5 Spatial variance spectra obtained for the 32^, 64^ and 128^ experiments 
(from top to bottom), normalized by a^z, where is the area below the spectral 
curve, calculated for a height z = 0.3h. The I, S and F chemical cases are rep- 
resented by solid, dashed and dash-dotted lines, respectively. The slope of -5/3 is 
also shown. 
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reacting scalars at those scales. As a consequence, in Fig. Elthe spectra cor- 
responding to the fast chemical regime show substantially smaller variance 
at the largest scales (small ki) for all the resolutions. In order to understand 
if that effect will extend to the subgrid and/or smallest resolved scales, one 
needs to consider, in addition, the value of the subgrid Damkohler number 
Dttsgs- If Ddsgs << 1 thcu thc subgrid and smallest resolved scales are 
not affected and, consequently, the scale-dependent dynamic model param- 
eters are not expected to depart from the inert scalar values, as shown in 
Fig. 21 On the other hand, if Daggs > 0(1), then the effect of the chemistry 
extends also to the subgrid and/or smallest resolved scales, which become 
less isotropic. It is important to note that the tuning-free dynamic model 
adjusts to the less isotropic behavior of those scales through a reduction of 
the model parameters as shown in Fig. 01 

4 Summary 

A scale-dependent dynamic model for the subgrid-scale fluxes of react- 
ing scalars is presented and used in large-eddy simulations of atmospheric 
boundary layer reacting flows. This model uses a dynamic procedure to cal- 
culate the lumped coefficient {ScJ^^ ^pC-s) ™ eddy-diffusion subgrid-scale 
model as a function of the dynamics of the resolved scales. Consequently, 
the proposed dynamic model does not require any parameter specification 
or a priori tuning. 

The scale-dependent dynamic model is implemented in simulations of a 
neutral ABL with a constant and uniform surface flux of a first-order decay- 
ing scalar. A set of spatial resolutions and reaction rates are used in order 
to study the ability of thc dynamic model to adjust the model coefficient 
to changes in scale (resolution) and/or chemical regime. The model is able 
to account for the scale dependence of the model coefficient in a self consis- 
tent way. Scale dependence increases with decreasing height-to-scale ratio 
due to the increased flow anisotropy at the fllter and/or test filter scales. 
In addition, we found that when the time scale associated with the subgrid 
turbulence is on the same order of magnitude that the time scale associated 
with the chemistry, i.e. Daggs ~ 0(1), an additional dependence towards the 
chemical reactivity can be expected. By calculating explicitly the lumped 
model coefficients for the reacting scalars, there is no need for any assump- 
tion on the combined effects of resolution and reaction rate on the efficiency 
of the subgrid scale turbulence to transport scalars. We showed that the 
model is able to capture the chemical regime dependence of the reactant 
and, as a result, simulations yield resolution-independent total turbulent 
reactant fluxes. 
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